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Abstract 

The abundance of gut microbiota can be viewed as a quantitative trait, wliicli is affected by the genetics and environment 
of the host. To quantify the effects of host genetics, we calculated the heritability of abundance of specific microorganisms 
and genetic correlations among them in the gut microbiota of two lines of chickens maintained under the same husbandry 
and dietary regimes. The lines, which originated from a common founder population, had undergone >50 generations of 
selection for high (HW) or low (LW) 56-day body weight and now differ by more than 10-fold in body weight at selection 
age. We identified families of Paenibadllaceae, Streptococcaceae, Helicobacteraceae, and Burkholderiaceae that had moderate 
heritabilities. Although there were no obvious phenotypic correlations among gut microbiota, significant genetic 
correlations were observed. Moreover, the effects were modified by genetic selection for body weight, which altered the 
quantitative genetic background of the host. Heritabilities for Bacillaceae, Flavobacteriaceae, Helicobacteraceae, 
Comamonadaceae, Enterococcaceae, and Streptococcaceae were moderate in LW line and little to zero in the HW line. 
These results suggest that loci associated with these microbiota families, while exhibiting genetic variation in LW, have been 
fixed in HW line. Also, long term selection for body weight has altered the genetic correlations among gut microbiota. No 
microbiota families had significant heritabilities in both the LW and HW lines suggesting that the presence and/or absence 
of a particular microbiota family either has a strong growth promoting or inhibiting effect, but not both. These results 
demonstrate that the quantitative genetics of the host have considerable influence on the gut microbiota. 
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Introduction 

High levels of diversity, community structure, and composition 
of gut microbiota are strongly associated with host species and are 
very stable and consistent within host species [1,2], reflecting 
selection of microorganisms by their host and environmental 
factors. Within host species, factors such as diet, litter, and 
maternal effects [3-6] , as well as single gene differences, affect the 
population structure of gut microbiota [7-10]. Previously, we 
reported that the quantitative genotype of the host influenced gut 
microbiota composition in adults from two lines of chickens fed a 
common diet and maintained under the same husbandry [11]. 
This finding suggests that the gut microbiome can be viewed as a 
complex, polygenic trait with a synergy of host genetic and 
environmental factors that shape and account for variability. UntU 
now, the role of host genetics on shaping this vital 'microbial 
organ' was not clear, and to our knowledge, lacking quantification 
of contributions from the polygenic background of the host. 



The genotype of the host may affect its microbiota composition 
directly through secretions into the gut, influences on gut motility 
and modification of epithelial cell surfaces, or indirectiy through 
food and lifestyle preferences. Detecting their role(s) wiU require 
well controlled effects other than those of the host genotype. Thus, 
choosing a model organism maintained in essentially an identical 
environment with few maternal effects should enhance our 
understanding of host genotype effects on gut microbiota. Here, 
we choose for our model two lines of chickens that had undergone 
>50 generations of bidirectional selection for 56-day high (HW) or 
low (LW) body weight (Figure 1). The lines originated from a 
common founder population and have complete pedigrees 
[12,13]. Throughout all generations, they have been maintained 
at the same location and reared on the same diets. Selection has 
resulted in more than a 10-fold difference between them for body 
weight at selection age. They are segregating populations with 
maximum inbreeding coefficients (F) being 0.53 and 0.61, with a 
mean of 0.26 (SD 0.15) and 0.30 (SD 0.17) in LW and HW lines. 



PLOS ONE I www.plosone.org 



1 



March 2014 | Volume 9 | Issue 3 | e89862 



Heritability of Gut Microbiota 



Bac" 



iWales 




B 



Bacillales 8% 



-n 




1 


Bacteria 


o 
c 

(0 


LW 



1 more 




'es 



iBacteroidaceae 
iMicrobacteriaceae 
iFlavobacterlaceae 
iRhodobacteraceae 
.L Qtnamonadaceae 
ampylobacteraceae 
LdTtobacillaceae 
ICorynebacteriaceae 
lEubacterlaceae 
iFusobacteriaceae 
iHelicobacteraceae 
[Micrococcaceae 
jrphyromonadaceae 
dsteurellaceae 
' arillaceae 

Viethanocorpusculaceae 
iDemococcaceae 
iBifidobacterlaceae 
ik enellaceae 
tieptococcaceae 
lAcidaminococcaceae 
Enterococcaceae 

?tococcaceae 
lEnterobacteriaceae 
lAcetobacteraceae 
l56weightlog 
IXanthomonadaceae 
ICellulomonadaceae 
iDesulfovibrionaceae 
iNocardioidaceae 
lAlcaligenaceae 
IRhlzobiaceae 
iMoraxellaceae 
ISphingobacteriaceae 
iBrevibacteriaceae 
lActlnomycetaceae 
ICarnobacterlaceae 
iRuminococcaceae 
IPeptostreptococcaceae 
lErysipelotrlchaceae 
ICaulobacteraceae 
iLachnospiraceae 
IPropionibacteriaceae 
IPeptococcaceae 
jPseudomonadaceae 
IVeillonellaceae 
IPIanococcaceae 
iBurkholderiaceae 
IPaenlbacillaceae 
IStaphylococcaceae 
Iciostridiaceae 




PLOS ONE I www.plosone.org 



2 



March 2014 | Volume 9 | Issue 3 | e89862 



Heritability of Gut Microbiota 



Figure 1 . High (HW) and low (LW) lines of chicken. A. The distribution of the gut nnicrobiomes composite for HW and LW lines. B. Some heritabilities 
of taxonomy in these lines. C. Picture of high (HW) and low (LW) line females. 
doi:1 0.1 371 /journal.pone.0089862.g001 



respectively [13]. QTL mapping revealed 13 loci afTecting growth 
in these two lines, however, each locus explained only a small 
additive effect for this large phenotypic difference [14]. This 
moderate F, along with QTL results, suggest that there is genetic 
diversity within each line, and the magnitude of the body weight 
difference between them at selection age is because of an 
accumulation of quantitative genes, each with very small effects 
rather than a single gene mutation. Not only do these lines provide 
an ideal model for studying quantitative genetic effects of the host 
on the microbiome of the gut, they also help us characterize how 
selection pressures on the host alter its genetic impact on the gut 
microbiota. Aided by next generation sequencing technology', we 
investigated the population structure of the gut microbiota in 
adults of these two lines. We treated the abundance of 
microorganisms of gut microbiota as a quantitative trait of the 
host by calculating the heritabUities of abundance of specific 
microorganisms and quantifying diversit)' and composition of 
community structure of the gut microbiome based on the host's 
quantitative genetic background. Genetic correlations were 
calculated among different microorganisms as they allow for 
quantifying the contribution from the host's genetic background to 
the interaction among the different microorganisms. 

Results and Discussion 

At least 22,000 sequence reads, which passed our quality 
criteria, were generated for each sample. The average sequence 
reads per sample was 42,358. Afti;r taxonomic classification, 50 
families, which were present in at least 65 samples, were 
considered as common and their normalized abundance counts 
[1 1] were used for further analysis. 

Selection pressure alters additive effects from the host's 
genetic background for the gut microbiome 

The abundance of microorganisms of gut microbiota, as a 
quantitative trait of the host, is affected by the additive effects from 
host genotype + environmental factors + random error. The 
heritability of a quantitative trait estimates the fraction of the 
phenotypic variation attributed to genetic variation. Here, we 
calculated the heritability of abundance of each microorganism to 
estimate the additive effects of the host's quantitative genetic 
background (Table 1). The heritability of 56-day body weight for 
the pooled group (combine the two lines and treat line and sex as 
fixed effects) was 0.326. This vEilue was consistent with estimates 
for this trait in other populations, and suggests that the moderate 
heritability of body weight for these populations is typical for that 
in chickens. Four families of microbiota [Paenibacillaceae, Strepto- 
coccaceae, Helicohackraceae, and Burkholderiaceae) had moderate 
heritabilities (Table 1). 

Although the two lines of chickens reported here originated 
from the same founder population, after >50 generations of long 
term pedigree selection there were changes in gene frequencies 
that could have altered their quantitative genetic background [14]. 
Their phenotypes (Figure 1) and physiological functions are quite 
different [13,15]. To determine if the changing quantitative 
genetic background of the host affected the additive effects for gut 
microbiota, we calculated heritabilities within each line, and the 
results were quite different (Table 1). Seventeen families had 
heritabilities greater than 0.20 in the LW line, 6 in the HW line, 
with the overlap between lines of only 2 families. Heritabilities of 



Bacillaceae, Flavobacteriaceae, Helicobacteraceae, Comamonadaceae, Enter- 
ococcaceae, and Streptococcaceae were moderate (p<0.05) in the LW 
hne and littie to zero in HW line. FamiUes of Aerococcaceae, 
Enterococcaceae, Lactobacillaceae, and Streptococcaceae belong to Lacto- 
bacillales, the order of most probiotics. Their heritabilities were 
moderate in LW and little to zero in HW. After >50 generations 
of selection for high or low body weight under a common diet, one 
may speculate that the quantitative genetic background is quite 
different in the two lines and genes which enhance growth in line 
HW have afforded sufficient absorption of nutrients so that over 
time there was less selection pressure for gut microbiota in HW. 
This suggests that for the LW line, there was greater need for these 
microorganisms, hence, greater effects from the host. This 
reasoning also supports the thesis that abundance of some 
microorganisms is determined by the genetic background of the 
host with the effect of each gene being small. 

Genetic correlations to quantify the association of host 
genetic background to microorganism interactions 

The genetic correlation is the proportion of variance (i.e. 
covariance) that two traits share due to common genetics. That the 
heritability of some microorganisms is moderate implies that they 
are influenced by the genotype of the host. Accordingly, we 
calculated genetic and phenotypic correlations among microor- 
ganisms in the gut microbiome whose heritabilities were greater 
than 0.20 (Table 2). Although there were no obvious phenotypic 
correlations among gut microorganisms, some genetic correlations 
were high, such as Streptococcaceae and Propionibacteriaceae (0.904**) 
and Paenibacillaceae and Bacillaceae (0.804**), while others such as 
Streptococcaceae and Comamonadaceae (—0.524**) and Helicobacteraceae 
and Comamonadaceae (0.480**) were moderate. We also calculated 
phenotypic and genetic correlations among gut microorganisms in 
die HW (Table SI) and LW (Table S2) lines, respectively. The 
phenotypic correlations between Burkholderiaceae and Paenibacillaceae 
were small in both lines (Table S2, Table S3), however, the genetic 
correlations were moderate (—0.442) in LW and high (—0.711**) 
in HW chickens. These results imply further the importance of the 
genetic background of the host on the interactions among 
microbiota. 

Correlations of 56-day body weight per se with gut 
microbiota 

Because 56-day body weight is a quantitative trait, genetic and 
phenotypic correlations between it and microorganisms whose 
heritability was greater than 0.20 were calculated (Table 2). 
Although no phenotypic correlation between 56-day body weight 
and microorganisms was significant, the genetic correlation with 
Deinococcaceae (0.513*) was positive while those for Paenibacillaceae 
(-0.449*), Bacillaceae (-0.433*), Enterococcaceae (-0.627**) and 
Acetobaeteraeeae (—0.524*) were negative (Table 2). To determine if 
these relationships were consistent after long term selection, we 
also calculated the phenotypic and genetic correlations between 
56-day body weight and gut microorganism in the HW (Table SI) 
and LW (Table S2) lines, respectively. Both Burkholderiaceae and 
Paenibacillaceae had moderate heritabilities in both hnes, however, 
their interaction with 56-day body weight was reversed. For 
Burkholderiaceae, the genetic correlation was —0.722 in HW, and 
0.583 in LW. For Paenibacillaceae, the genetic correlation was 0.543 
in FIW and —0.747 in LW, implying further that no major single 
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Table 1. Heritabilities by line and pooled. 
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Pooled 


Low Weight 


High Weight 


Acetobacteraceae 


0.196 


0.216 


0.181 


Methanocorpusculaceae 


0.080 


0.210 


0.000 



*: p<0.05, 
**: p<0.01. 

doi:1 0.1 371 /journal.pone.0089862.t002 



gene was responsible for host-microorganism interactions. Rather, 
they are influenced by several genes each witli small eflects and 
can be bidirectionally dependent on the host. 

General 0043omments 

Heritabilities of Bacillaceae, Flavobacteriaceae, Helicobacteraceae, 
Comamonadaceae, Enterococcaceae, and Streptococcaceae were moderate 
(p<0.05) in the LW line and htde to zero in the HW line. These 
results suggest that genes which contribute to the abundance of 
these microorganisms families may be frxed in the HW but not in 
the LW line. It also can be deduced that the presence and/or 
absence of these microbe families are more important in the HW 
line because the heritability is 0 implying no genetic variation. 
Because the heritability of these families are greater than 0 in one 
of the lines (LW) then genetic variation contributing to these 
microbiota families was likely present in the founding population. 
There is a novel deletion that is fixed in the HW line and occurs at 
a low frequency in the LW line, which contributes to the 
diflerences in growth [16] and feed efficiency [Siegel, unpub- 
lished]. Combined with our results, it suggests that some genes 
located in the deletion region may be associated with microbiota 
families. 

None of the microbiota families that we measured had 
significant heritabilities in both the LW and HW lines, suggesting 
that the presence and/ or absence of a particular microbiota family 
can either have a strong growth promoting or inhibiting effect, but 
not both. 

In order to avoid reproduction problems and metabohc issues 
caused by obesity, the HW line chickens used in this experiment 
were feed restricted starting at 56 days, which was selection age. 
This practice initiated in generation SI 8, was necessary to prevent 
obesity from overconsumption from ad libitum feeding. Because 
the objective of this experiment was to address if the genetic 
background of the host influenced the gut microbiome commu- 
nity, we collected fecal samples at 245 days of age when the host 
was mature. Studies show that dietary composition and caloric 
intake appear to swiftly regulate intestinal microbial composition 
and function [17,18]. Therefore, fecal samples were taken 189 
days after feed restriction was started, a period that should have 
allowed the chickens to sufficiendy adjust to issues associated with 
feed restriction. 

Our results are consistent with the thesis that the genetic 
background of the host influences the diversity and composition of 
gut microbiome. Involved are numerous genes, each having a 
small additive effect. Moreover, the quantitative genetic expression 
of the host win change over time due to selection, resulting in 
different families of microorganisms. Although gut microbiota 
could influence physiological functions of the host, and vice versa, 
its accumulation in the gut is greatly influenced by the host's 
genotype. 
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Conclusions 

Families of Paenibacillaceae, Streptococcaceae, Helicobacteraceae, and 
BurkhoUeriaceae had moderate heritabilities, which provide evidence 
that the abundance of gut microbiota can be viewed as a 
quantitative trait and evidence of an association between host 
genome and gut microbiota. Moreover, these associations have 
been modified by genetic selection for body weight, which altered 
gene frequencies of the host resulting in a change in its quantitative 
genetic background. The outcome of the divergent selection is that 
heritabilities for BaciUaceae, FhwobacUriaceae, Helicobacteraceae, Coma- 
monadaceae, Enterococcaceae, and Streptococcaceae were moderate and 
litde to zero in the LW and HW lines, respectively. 

Although there were no obvious phenotypic correlations among 
gut microbiota, there were significant genetic correlations, imply- 
ing that the long term selection for body weight also altered the 
genetic correlations among gut microbiota. Therefore some of 
interactions among gut microbiota are mediated by the host. 
These mediation are through host genes, and changing these gene 
frequencies will also modify interactions among gut microbiota. 
Thus, the quantitative genetic background of the host can have 
considerable influence on the gut microbiota. 

Materials and Metliods 

Animals and sample collection 

Protocols used for this experiment were approved by the 
Institutional Animal Care and Use Committee at Virginia Tech. 
Fresh fecal samples from 258 adult chickens, which consisted of 88 
LW females, 38 LW males, 90 HW females, and 42 HW males, 
were collected. The fecal samples were obtained from generations 
54 and 55 of the HW and LW lines when the chickens were 245 
days of age. These chickens had been under similar husbandry 
conditions and were fed a corn-soybean, non-pathogen free diet in 
mash form. Details are provided in the supplement of breeder diet 
ingredients. Cages had sloping wire floors with papers beneath 
them to collect feces. 

DNA extraction, PCR amplification of 16S rRNA, amplicon 
sequence and sequence data processing 

Microbial genome DNA was extracted from fecal samples using 
QIAamp DNA stool mini kit (QIAGEN, cat#51504) following the 
manufacturer's recommendation. According to our previous 
comparisons, results from V3, V4, V1-V3, and V4— V6 hyper- 
variable regions of 16S rRNA [1 1], V4 of 16S rRNA were PCR 
amplified from microbial genome DNA harvested from fecal 
samples using barcoded fusion primers (forward primers: 
5'AYTGGGYDTAAAGNG 3', reverse primers: 5'TACNVG- 
GGTATCTAATCC 3') and were used for the remainder of our 
study. The PCR condition and PCR product purification followed 
our previous publication [1 1] . Barcoded V4 amplicons were 
sequenced using the pair-end method by lUumina Miseq with a 6 
cycle index read. Sequences with an average phred score lower 
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